Bumpy Black Holes in Alternative Theories of Gravity 
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We generalize the bumpy black hole framework to allow for alternative theory deformations. We 
construct two model-independent parametric deviations from the Kerr metric: one built from a 
generalization of the quasi-Kerr and bumpy metrics and one built directly from perturbations of the 
Kerr spacetime in Lewis-Papapetrou form. We find the conditions that these "bumps" must satisfy 
for there to exist an approximate second-order Killing tensor so that the perturbed spacetime still 
possesses three constants of the motion (a deformed energy, angular momentum and Carter constant) 
and the geodesic equations can be written in first-order form. We map these parameterized metrics to 
each other via a diffeomorphism and to known analytical black hole solutions in alternative theories of 
gravity. The parameterized metrics presented here serve as frameworks for the systematic calculation 
of extreme mass-ratio inspiral waveforms in parameterized non-GR theories and the investigation 
of the accuracy to which space-borne gravitational wave detectors can constrain such deviations. 

PACS numbers: 04.30.-w,04.50.Kd,04.70.-s 



I. INTRODUCTION 

Gravitational waves (GWs) will be powerful tools for 
learning about source astrophysics and testing strong 
field gravity [iLM- The Laser Interferometer Space An- 
tenna (LISA) 0t9, for example, is expected to be sen- 
sitive to roughly a few radians or better of GW phase 
in a one-year observation. Such levels of precision are 
achieved through matched filtering, where the data is 
cross-correlated with a family of waveform models (see 
eg. Q)- If a signal is present in the data, this cross- 
correlation acts as a filter that selects the member of the 
waveform family that most closely resembles the signal. 

Extreme mass-ratio inspirals (EMRIs) are ideal astro- 
physical sources to perform precision GW tests of strong 
field gravity with LISA [H, % . These sources consist of 
a small compact object (SCO) that spirals into a super- 
massive black hole (SMBH) in a generic orbit, producing 
millions of radians in GW phase inside LISA's sensitiv- 
ity band. These waves carry detailed information about 
the spacetime geometry in which the SCO moves, thus 
serving as a probe to test general relativity (GR). 

Tests of strong field gravity cannot rely on waveform 
families that assume GR is correct a priori. Instead, such 
tests must employ more generic waveforms that allow for 
GR deviations. Recently, Yunes and Pretorius pro- 
posed the parameterized post-Einsteinian (ppE) frame- 
work, in which analytic waveforms that represent com- 
parable mass-ratio coalescences in GR are parametrically 
deformed. For specific values of the ppE parameters, one 
recovers GR waveforms; otherwise they describe wave- 
forms in non-GR theories. The ppE scheme has been 
shown to be sufficiently flexible to map all known al- 
ternative theory predictions for comparable mass-ratio 
coalescences 

PpE deformations of EMRI waveforms have not yet 
been constructed because of the complexities associated 
with the computation of such waveforms. EMRIs are 
not amenable to post-Newtonian (PN) approximation 



schemes EMI, 

from which many GR templates arc to- 
day analytically constructed. Instead, to build EMRI 
waveforms one must solve geodesic-like equations, en- 
hanced with a radiation-reaction force that induces an 
inspiral. The solution of such differential equations can 
only be found numerically (see eg. [l^ for a recent review 
of the self- force problem). 

One can parametrically deform EMRI waveforms by 
introducing non-GR deviations in the numerical scheme 
used to build such waveforms. Two main ingredients can 
then be modified: the conservative sector, which controls 
the shape of non-radiative orbits, and the dissipative sec- 
tor, which controls the rate of inspiral and the GW gen- 
eration mechanism. The conservative sector, to leading 
order in the mass-ratio, depends only on the background 
spacetime metric, the Kerr metric, on which the SCO 
moves like a test-particle. 

The goal of this work is to find a parametric deforma- 
tion of the Kerr metric that allows for non-GR deviations 
while retaining a smooth Kerr limit, i.e. as the deforma- 
tion parameters go to zero, the deformed Kerr metric 
reduces smoothly and exactly to Kerr. This limit then 
guarantees that the main properties of the Kerr back- 
ground survive to the deformation, such as the existence 
of an event horizon, an ergosphere, and, in particular, 
three constants of the motion, such that the geodesic 
equations can be separated into first-order form. 

Although motion in alternative theories of gravity need 
not be geodesic, we will here assume that the geodesic 
equations hold to leading order in the mass ratio. This 
then allows us to focus only on metric deformations that 
correct the conservative sector of EMRI waveforms. This 
assumption is justified for theories that derive from cer- 
tain diffeomorphism-covariant Lagrangians [l6j , as is the 
case, for exarnple, in dynamical Chern-Simons (CS) mod- 
ified gravity [131, i'^ Einstein-Dilaton-Gauss-Bonnet the- 
ory [l8l - |22j . and in dynamical quadratic gravity [23l - [25| . 
We are not aware of any alternative theory where this 
is not the case to leading order in the mass ratio (and 
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neglecting spins), although it is not hard to imagine that 
one may exist. 

A framework already exists to parametrically deform 
the metric tensor through the construction of so-called 
bumpy spacetimes j26l428{ . Such metrics are deforma- 
tions of the Kerr metric, where the bumps are required 
to satisfy the Einstein equations. Because of this last 
condition, bumpy black holes (BHs) are not sufficiently 
generic to represent BHs in certain alternative theories of 
gravity (e.g. solutions that are not Ricci flat). A better 
interpretation is to think of these bumps as representing 
exterior matter distributions. The bumpy BH formalism 
then allows tests of whether compact objects are truly 
described by the vacuum Kerr metric or by some more 
general tensor with external matter sources, assuming 
GR still holds. 

We here propose two generalizations of the bumpy BH 
formalism to allow for metric deformations that can rep- 
resent vacuum BH solutions in alternative theories. The 
first approach, the generalized bumpy Kerr (BK) scheme, 
takes the standard bumpy metric and relaxes the require- 
ment that the bumps satisfy the Einstein equations. The 
second approach, the generalized deformed Kerr (DK) 
metric, perturbs the most general stationary, axisymmet- 
ric metric in Lewis-Papapctrou form and transforms it to 
Boycr-Lindquist coordinates without assuming the defor- 
mations satisfy the Einstein equations. In both cases, the 
metric is exactly Kerr to zeroth order in the deforma- 
tion parameter, but it deviates from Kerr at first order 
through a set of arbitrary functions. 

The metric deformation is then restricted by requir- 
ing that the full metric still possesses three constants of 
the motion: a conserved energy (associated with time 
invariance), a conserved angular momentum component 
(associated with azimuthal rotational invariance), and an 
approximate second-order Killing tensor that leads to a 
conserved Carter constant. These conditions restrict the 
set of arbitrary functions that parameterize the metric 
deformations to have a certain functional form. It is these 
conditions that select the appropriate "bumps" instead 
of the imposition of the Einstein equations. 

The restrictions imposed above are not strictly nec- 
essary, as there is no guarantee that BHs in alternative 
theories of gravity will continue to have three conserved 
quantities. However, all BH solutions in alternative the- 
ories of gravity known to date and that are not patho- 
logical (i.e. they are stationary, axisymmetric, asymptot- 
ically fiat, and contain no spacetime regions with closed 
time-like curves outside the horizon) possess three con- 
stants of the motion [2^ . Examples include the slowly- 
rotating solution found in dynamical CS gravity (sOj and 
the spherically symmetric solution found in dynamical 
quadratic gravity (25| . 

Much effort has gone into finding spacetimes with an 
exact second-order (or higher-order) generalized "Carter 
constant," thus allowing for the separation of the equa- 
tions of motion [3ll - l35| . That work demonstrates that, 
for a broad class of spacetimes, such separation can 



be done provided that the Carter constant is quartic 
in the orbit's 4-momentum (i.e., the constant is C = 
£.ai3-ySP"p^p''p^ , where (.a/SjS is a 4th-rank Killing tensor 
and p" is the 4-momentum). In this analysis, we show 
that one can in fact find an approximate Carter constant 
that is quadratic in the 4-momentum, C — £,app"p^ , for 
many relevant spacetimes, provided they differ from Kerr 
only perturbatively. 

Finally, we gain insight on the different proposed pa- 
rameterizations by studying certain key limits. First, we 
show that the BK and DK deformations are related to 
each other by a gauge transformation. Second, we show 
that both the BK and DK metrics can be exactly mapped 
to specific non-GR BH metrics, the dynamical CS gravity 
one and the dynamical quadratic one. Third, we study 
the structure that the deformations must take when only 
deforming frame-dragging or the location of the event 
horizon. Fourth, we separate the geodesic equations into 
first-order form in both the BK and DK metrics, and cal- 
culate the modified Kepler law. This modification cor- 
rects the dissipative dynamics through the conversion of 
radial quantities to frequency space. 

The parameterized metrics (BK and DK) proposed in 
this paper lay the foundations for a systematic construc- 
tion and study of ppE EMRI waveforms. With these 
metrics and their associated separated equations of mo- 
tion, one can now study modified SCO trajectories and 
see how these impact the GW observable. Numerical im- 
plementation and a detailed data analysis study will be 
presented in a forthcoming publication. 

This paper is organized as follows. Sec HIl reviews BHs 
in GR and derives the first-order form of the geodesic 
equations. Sec. IIIII introduces the BK formalism, which 
generalizes the bumpy BH formalism. Sec. IIVI presents 
the DK parameterized metric. Sec. |V] compares the pa- 
rametcrizations to each other and to alternative theory 
predictions. Sec. IVII discusses some general properties of 
the parameterizations and points to future work. 

In this paper we primarily follow the notation of Mis- 
ner, Thorne and Wheeler [3Q|. Greek letters stand for 
spacetime indices, while Latin letters in the middle of the 
alphabet (i, j, k, . . .) stand for spatial indices only. Paren- 
thesis and brackets in index lists stand for symmetriza- 
tion and antisymmetrization respectively, i.e. A(^g_i,-^ = 
{Aab + Aba)/2 and A[ab] = {Aab - Aba)/2. Background 
quantities are denoted with an overhead bar, such that 
the full metric g^i^ = g^i^ -I- eh^^i, can be decomposed into 
a background (jfn, plus a small deformation h^^^, where 
e ^ 1 is a book-keeping parameter. The background 
spacetime is always taken to be the Kerr metric. We use 
geometric units with G = c = I. 



II. BLACK HOLES AND TEST-PARTICLE 
MOTION IN GENERAL RELATIVITY 

In this section, we review some basic material on the 
Kerr BH solution in GR and the motion of test-particles 
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(see e.g. [S^ for more details). The most general station- 
ary, axisymmetric BH solution in GR is the Kerr metric, 
which in Boyer-Lindquist coordinates (t, r, 6, (/>) takes the 
form 
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where we have introduced the quantities 



+ cos^ Q . 



A = r^f + a\ 
S = (r^ +0^)2 - a^Asin^e* 
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(1) 

(2) 
(3) 
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The overhead bars are to remind us that the Kerr metric 
will be taken as our background spacetime. This back- 
ground describes a BH with mass M and spin angular 
momentum 5* ~ (0,0, A/a), where a is the Kerr spin 
parameter with units of length. 



A. Equations of Motion 

With this background, we can now compute the equa- 
tions of motion for a test particle of mass /z moving on 
a background worldline i'' = ^''(A), where A is an affine 
parameter. One way to derive these equations is via the 
action for a non-spinning test-particle (see e.g. [37| ) 



5. 



(5) 



where = dz'^/dX is the tangent to and g is the 
determinant of the background metric. 

The contribution of this action to the field equations 
can be obtained by varying it with respect to the metric 
tensor. Doing so, we obtain the stress-energy for the 
test-particle: 

T^ii^n = / ^u"S^5W[a, - z{r)] , (6) 
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where 
dA 



r is proper 



\l~9o.fs{z)i 



time. 



related to A via dr ~ 
= dz^/dr is the parti- 
cle's 4- velocity, normalized via g^^u^W^ ~ —1 and 

6^^) is the four-dimensional Dirac density, defined via 
/d4xV=^5W(a;) = 1. 

The divergence of this stress-energy tensor must vanish 
in GR, which implies that test-particles follow geodesies: 



D dz" 
dr dr 



0. 



(7) 



where D/dr is a covariant derivative. The divergence 
of the stress-energy tensor vanishes due to local energy- 
momentum conservation and the equivalence princi- 
ple [1^. We have assumed here that the particle is non- 
spinning, so that geodesic motion is simply described by 



Eq. (O; otherwise other terms would arise in the action 
that would lead to spin-dependent modifications. 



B. Constants of the Motion and Separability 

The geodesic equations as written in Eq. ([7]) are of 
second-order form, but they can be simplified to first- 
order form if there exist at least three constants of 
the motion plus a normalization condition on the four- 
momentum. If the metric is sufficiently symmetric, then 
Killing vectors ^" and Killing tensors ^"'^ exist that sat- 
isfy the vectorial and tensorial Killing equations 



V(a^/3) = 0, 



(8) 
(9) 



The Kerr metric possesses two Killing vectors and 
one Killing tensor. The vectors, = [1,0,0,0] and 
= [0, 0, 0, 1], are associated with the stationarity and 
axisymmetry of the metric respectively. In addition, the 
Kerr metric also possesses the Killing tensor: 



(afi = Akf^Jp) +r g^i3, 
where fc" and P are principal null directions: 
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The contraction of these Killing vectors and tensors 
allows us to define constants of the motion, namely the 
energy E, z-component of the angular momentum L, and 
the Carter constant C: 



E 
L 
C 



(12) 
(13) 
(14) 



It is conventional to define another version of the Carter 
constant, namely 

\2 



Q = C-{L-aEy 



(15) 



With these three constants of the motion [either 
(E, L, C) or [E, L, Q)] and the momentum condition 
p"Pa ~ ^M^, or simply Uqu" = — 1, we can separate the 
geodesic equations and write them in first-order form: 
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where overhead dots stand for differentiation with respect 
to proper time and P = E{r^ + a^) — aL. Notice that not 
only have the equations been written in first-order form, 
but they have also been separated. 

Once the geodesic equations have been derived, one can 
compute Kepler's third law, relating radial separations to 
orbital frequencies for an object moving in an equatorial, 
circular orbit. For such an orbit, the equatorial, circular 
orbit conditions f = and df/dr = lead to [38| 

In the far field limit, M/r ^ 1, this equation reduces to 
the famihar Kepler law: fl ^ (M/r^)^/^. 



III. BLACK HOLES AND TEST-PARTICLE 
MOTION IN ALTERNATIVE THEORIES: 
GENERALIZED BUMPY FORMALISM 

We here introduce the standard bumpy BH formal- 
ism (26l - [28j and generalize it to the BK framework. The 



bumps are constrained by requiring that an approximate, 
second-order Killing tensor exists. We then rewrite the 
geodesic equations in first-order form. 

A. From Standard to Generalized Bumpy Black 
Holes 

Let us first review the basic concepts associated with 
the bumpy BH formalism [26l - l28| . This framework was 
initially introduced by Collins and Hughes [l^ to model 
deformations of the Schwarzschild metric and expanded 
by Vigeland and Hughes [23] to describe deformations on 
a Kerr background. In its standard formulation, gener- 
alized to spinning BHs [^^j, the metric is perturbatively 
expanded via 

gf^iy =gf,i^ + <^ h'jfj' , (21) 

where e is a book-keeping parameter that reminds us that 
\h^^J^\/\gfj,u\ ^ 1. The background metric g^i, is the Kerr 
solution of Eq. ([T|), while in the standard bumpy formal- 
ism, the metric functions h^^^ are parameterized via 
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(22) 



where {ipl^^ , Ji^^ , crf^^) are functions of {r,9) only that 
represent the "standard bumps." 

A few comments are due at this point. First, notice 
that this metric contains only three arbitrary functions 
{ipf^'^,ji^'^,ai°'^), instead of four, as one would expect 
from the most general stationary and axisymmetric met- 
ric. This is because of the specific way Eq. ((22|) is de- 
rived (see below), which assumes a Ricci-flat metric in 
the a = limit. Second, note that many metric com- 
ponents are non- vanishing: apart from the usual (t,t), 
(t, 0), (r, r), (9,0) and (0,0) metric components that 
are non-zero in the Kerr metric, Eq. also has non- 
vanishing (i, r) and (r, 0) components. 

The derivation of this metric is as follows. One begins 
with the most general stationary and spherically symmet- 
ric metric in Lewis-Pappapctrou form that is Ricci-flat. 
One then perturbs the two arbitrary functions in this 
metric and maps it to Schwarschild coordinates. Through 



the Newman- Janis procedure, one then performs a com- 
plex rotation of the tetrad associated with this metric 
to transform it into a deformed Kerr spacetime. At no 
stage in this procedure is one guaranteed that the result- 
ing spacetime will still possess a Carter constant or that 
it will remain vacuum. In fact, the metric constructed 
from Eq. docs not have a Carter constant, nor does 
it satisfy the vacuum Einstein equations, leading to a 
non-zero effective stress-energy tensor. 

Such features make the standard bumpy formalism not 
ideal for null-tests of GR. One would prefer to have a 
framework that is generic enough to allow for non-GR 
tests while still possessing a smooth GR limit, such that 
as the deformations are taken to zero, one recovers ex- 
actly the Kerr metric. This limit implies that the de- 
formed metric will retain many of the nice properties of 
the Kerr background, such as the existence of an event 
horizon, an crgosphere, and a Carter constant. 
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This can be achieved by gencrahzing the bumpy BH 
formahsm through the relaxation of the initial assump- 
tion that the metric be Ricci-flat prior to the Newman- 
Janis procedure. Inspired by [26"28!], we therefore pro- 
mote the non-vanishing components of the metric per- 
turbation, i.e., ihf,^,hf:.\hf^,K^,h^^^,hlf,h''^^), to ar- 
bitrary functions of (r, 9). These functions are restricted 
only by requiring that an approximate second-order 
Killing tensor still exists. 

This restriction is not strictly necessary, but it is ap- 
pealing on several fronts. First, the few analytic non-GR 
BH solutions that are known happen to have a Carter 
constant, at least perturbatively as an expansion in the 
non-Kerr deviation. Second, a Carter constant allows for 
the separation of the equations of motion into first-order 
form, which then renders the system easily integrable 
with already developed numerical techniques. 



B. Existence Conditions for the Carter Constant 

Let us now investigate what conditions must be en- 
forced on the metric perturbation so that a Killing 
tensor and its associated Carter constant exist, at 
least perturbatively to 0{e). Killing's equations can be 
used to infer some fairly general properties that such 
a tensor must have. First, this tensor must be non- 
vanishing in the same components as the metric pertur- 
bation. For the generalized bumpy metric, this means 
that (5te, ^re, Ce^) must all vanish. Furthermore, if we 
expand the Killing tensor as 



(23) 



we then find that S^ap must have the same parity as 
^ai3 for the full Killing tensor to have a definite par- 
ity. For the generalized bumpy metric, this means that 
Crr, S.ee, ^00, 60) must be even under reflection: 9 -> 

e-n. 

These conditions imply that of the 10 independent de- 
grees of freedom in (5^q^, only 7 arc truly necessary. With 
this in mind, we parameterize the Killing tensor as in 
Eq. ((ini), namely 



(24) 



Notice that ^ap depends here on the full metric gap and 
on null vectors fc" and P that are not required to be 
the Kerr ones or the principal congruences of the full 



spacetime: fc" 7^ fc" and 7^ P. We decompose these 
vectors via 

fc" = fc" + e (5fc" , r = P + e Sr , (25) 

where fc" and P are given in Eq. pT|) . and the tensor 
into 



where ^q,^ is given in Eq. (1101) . while 
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(27) 



All lowering and raising of indices is carried out with the 
background metric. Although this particular ansatz only 
allows 6 independent degrees of freedom (as the vectors 
are assumed to be null), we will see it suffices to find a 
Carter constant. 

With this ansatz, the tensor Killing equation [Eq. ^] 
becomes 



(28) 



The system of equations one must solve is truly 
formidable and in fact overconstrained. Equation (j28p 
is a set of 20 partial differential equations, while the nor- 
malization conditions l"la = = k°'ka add two addi- 
tional algebraic equations. This means there are a total 
of 13 degrees of freedom (7 in and 6 in Sl°' and Sk" 
after imposing the normalization condition) but 20 par- 
tial differential equations to solve. 

In spite of these difficulties, we have solved this system 
of equations with Maple and the GrTcnsorll package [s^ 
and found that the perturbation to the null vectors must 
satisfy 



^6kl^ + 5,,6kl^,0,-Skl,, + 5, 



— — 6kl^ + 63,6kl^ + (54, 0, -^Sk^^ + 65 



(29) 



where di = Si{r) are arbitrary functions of r, generated 
upon solving the differential system. These functions are 
fully determined by the metric perturbation, which is 
given by 
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where 



a^cos^O, 81.2 = 6i.2(^) are functions 



of 9 and _Pf'^ are polynomials in r and cos0, which are 
given explicitly in Appendix [X] Notice that Sk'^^^ does 
not appear in the metric perturbation at all, so we are 
free to set it to zero, i.e. (Jfe^^ = 0. 

With this at hand, given some metric perturbation 
that satisfies Eq. (PU)) . one can construct the arbitrary 
functions Si and 0,; and thus build the null directions 
of the perturbed spacetime, such that a Killing tensor 
and a Carter constant exist perturbatively to 0{e). We 
have verified that the conditions described above satisfy 
the generic Killing tensor properties, described at the 
beginning of this subsection. 



1. Non-Rotating Limit 

Let us now investigate the non-rotating limit a — > 0. 
To do so, we expand all arbitrary functions via 



(31) 



where JFi is any of Si (r) or Qi (6) . Let us also expand the 
metric components in a similar fashion, i.e. 

hTpir. 0) = hlloir, 0) + a hll,{r, 6) + ©(a^) . (32) 

The following choices 



©1.0 = = 82, „ , ^1,0 = 



(52.0 = = ^5,0 , 52.1 = S^,i = - 



/ BK LBK 

r - 4M , 



2^3/2 "tt,0' 



1^4,0 — / Kr.O ' '^3,0 



LBK LBK 

"-rr.O , "tt,0 



2/2 



(33) 



where we recall that / = 1 — 2M/r, force the metric 
components to take the form 

hl^ = hfM, Kf = Kf^,{r), (34) 

where we have set all integration constants to zero. All 
other components of the metric vanish to this order. To 



obtain this result, it is crucial to use the slow-rotation 
expansion postulated above, as some of the conditions 
in Eq. ([50]) have seemingly divergent pieces. We see then 
that in the non-rotating limit, one can always choose free 
functions 5i and Qi such that the only two independent 
metric perturbations are h^^ and h^^. 

The Killing tensor is then given by Eq. ([M]), with the 
perturbed null vector components 
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^,//C„,0,O(a) 



where we have set ^fc^j^ = 0. Composing the Killing 
tensor, we find that all 0(e) terms vanish in the a — > 
limit. This result makes perfect sense, considering that 
the Schwarzschild Killing tensor, ^q,^ = r'^O.^p with 
Q,ap = diag(0, 0, r^, sin^ 6), is also a Killing tensor for 
the most generic static and spherically symmetric space- 
time with arbitrary it, t) and (r, r) components. 



Equations of Motion and Geodesies in 
First-Order Form 



The equations of motion for test-particles in alterna- 
tive theories of gravity need not be geodesic. But as 
we review in Appendix [B] they can be approximated as 
geodesies in the test-particle limit for a wide class of al- 
ternative theories, provided we neglect the spin of the 
small body. Here, we restrict attention to theories where 
the modified equations of motion remain geodesic, but 
of the generalized bumpy Kerr background constructed 
above. 

The geodesic equations can be rewritten in first- 
order form provided that three constants of the mo- 
tion exist. In Sec. IIIIBl we showed that provided 
the generalized bumpy functions satisfy certain condi- 
tions, then an approximate, second-order Killing ten- 
sor and a conserved Carter constant exist. Further- 
more, it is easy to see that the metric in Eqs. (PT|) 
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and (|22)) remains stationary and axisymmctric, since 
{hft,hf:^,hf^,h-^,h-lh-f,h--^) are functions of only 

ir,e). 

The constants of the motion then become 



E = E + eih^^tf^u" +gf,af'Su'') , 
L = L + e {h^,<P<'u'' + g^.cjy^Su-') , 

C = C + e(^5^f,^u''u'' + 2^^^u^^'6u•''> 
while the normalization condition becomes 



(36) 
(37) 
(38) 



(39) 



since by definition — 1 = gf^^u'^u'^. In these equations, 

= [t, f, 9, (j)] is the unperturbed, Kerr four- velocity, 
while Su'^ is a perturbation. 

Equations (p6|) - ((39l) can be decoupled once we make 
a choice for {E,L,Q). This choice affects whether the 
turning points in the orbit are kept the same as in GR, 
or whether the constants of the motion are kept the same. 
We here choose the latter for simplicity by setting E = E, 
L = L and C = C, which then implies that Su'^ must 
be such that all terms in parenthesis in Eqs. (|36l) - ((55)) 
vanish. Using this condition and Eq. we then find 

that 

pH = TK{r,e)+dT{r,0), 



RK{r) + SR{r, 6*), 

p^e^ = eK{0) + se{r,e), 



(40) 



where {Tk,Rk,Qki^k) are given by the right-hand 
sides of Eqs. (HH]), (dH), and ^ respectively. The 

perturbations to the potential functions {6T, 6R, 60, (5$) 
are given by 



6T{r,e) = 



A 

2aMr, 



sin^ ( 



T-taU 



(5$(r,( 

5R{r, ( 
(5e(r, ( 



A 

2aMr 



2Mr , 



A '^'"^ Asin^fl 
A[A{r,0)r^ + B{r,e)] , 

A{r,e) cos^e- B{r, 6) 



(41) 

(42) 

(43) 
(44) 



where the functions A{r,9) and B{r,9) are proportional 
to the perturbation: 



hntt 



A{r,e) ^ 2 



(45) 



(46) 



As before, we have dropped the superscript BK here. 
Interestingly, notice that the perturbation automatically 



couples the (r, 9) sector, so that the first-order equations 
are not necessarily separable. One can choose, however, 
the 7i functions in such a way so that the equations re- 
main separable, as will be shown elsewhere. 

The perturbations to the potential functions result in 
a modification of Kepler's law [Eq. pO|) ]: 



\n\ = 



Ml/2 (^i/2(,._3M) + 2aA/i/2) 
\n\ ^ ^ , '-5T{r) 



j.5/4 (^3/2 + aMi/2) 
(ri/2(r_3A/) + 2aA/i/2)i/' 



^5/4 (^3/2 _^ aMl/2) 



(47) 



where we recall that was already given in Eq. (pO|) . ST 
and are given by Eq. (|^T|) and (|35|), and we have used 
that E = E tmAL^L. In the far field hmit {M/r < 1), 
this becomes 



1 



-(5$ 



(48) 



where we have assumed that 5$ and 5T are both of order 
unity for simplicity. Given any particular metric pertur- 
bation, one can easily recalculate such a correction to 
Kepler's law from Eq. (|47|. Clearly, a modification of 
this type in the Kepler relation automatically modifies 
the dissipative dynamics when converting quantities that 
depend on the orbital frequency to radius and vice- versa. 



IV. BLACK HOLES AND TEST-PARTICLE 
MOTION IN ALTERNATIVE THEORIES: 
DEFORMED KERR FORMALISM 

In the previous section, we generalized the bumpy BH 
framework at the cost of introducing a large number of 
arbitrary functions, later constrained by the requirement 
of the existence of a Carter constant. In this section, we 
investigate a different parameterization (DK) that iso- 
lates the physically independent degrees of freedom from 
the start so as to minimize the introduction of arbitrary 
functions. 



A. Deformed Kerr Geometry 

Let us first consider the most general spacetime met- 
ric that one can construct for a stationary and axisym- 
mctric spacetime. In such a geometry, there will exist 
two Killing vectors, t° and that represent invariance 
under a time and azimuthal coordinate transformation. 
Because these Killing vectors are independent, they will 
commute satisfying t^a'Pb^ ctd] = = t^a4'b^ cipd]- Let us 
further assume the integrability condition 



f'RjH'' 



= (l)''Rah''cj)'^^ 



(49) 



The condition in Eq. (j^Hl) guarantees that the 2-planes 
orthogonal to the Killing vectors t° and 0" arc integrable. 
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Generic stationary and axisymmetric solutions to modi- 
fied field equations do not need to satisfy Eq. P^ . How- 
ever, all known analytic solutions in GR and in alterna- 
tive theories do happen to satisfy this condition. We will 
thus assume it holds here as well. 

Given these conditions, the most general stationary 
and axisymmetric line element can be written in Lewis- 
Papapctrou canonical form, using cylindrical coordinates 
{t,p, (t),z): 

ds^ = -V {dt - wd(f>f + V-^p^dcf^ + [dp^ + Kdz^) , 

(50) 

where fi,A) are arbitrary functions of (/c, 2;). In 

GR, by imposing Ricci-flatness, one can eliminate one of 
these functions via a coordinate transformation, and we 
can thus set A = 1. The Einstein equations then further 
restrict the form of these arbitrary functions [i^ 



V = 



{r+ +r-y- 4Af 2 + (r+ - r^f 



2aM M 



1 _ (?-+-r_)- 
4(M2-o2) 



A = 1, 

n = V— -° ^ + -,(51) 

4r+r_ 



where r± = y -|- (z ± V — a^) . One can map 
from this coordinate system to Boyer-Lindquist coordi- 
nates (t,r,9,(j}) via the transformation p = -y/A sin 9 and 
z ~ (r — Af) sin 6* to find the metric in Eq. ([T]). 

We will not impose Ricci-fiatness here, as we wish 
to model alternative theory deviations from Kerr. We 
thus work directly with Eq. ([SO)) , keeping (F, tu, f2,A) 
as arbitrary functions, but transforming this to Boycr- 
Lindiquist coordinates. Let us then re-parametrize 51^ = 
7, A = 7A, and q = Vw, and perturb the metric functions 



V = V + SV, q=^q + 6q, 
A = A + (5A, 7 = 7 + (57. 

The metric then becomes 

where the perturbation is given by 

-SV, hf^^Sq, 

(r - A/)2 sin^ t 



(52) 
(53) 



. DK 

'ht 



. DK 

"■re 

. DK 



= 5A( 



(57- 



(r - M) cos e sin {S-f - SX) , 
S\ sin^ e{r- Mf + (57 cos^ 6IA , 

{AaM r Sq 



2Mr 



sin^ 0A 



5V 



}, (54) 



and all other components vanish. 

As is clear from the above expressions, the most general 
perturbation to a stationary, axisymmetric metric yields 
five non-vanishing metric components that depend still 
only on four arbitrary functions {SV, Sq, S^, SX) of two 
coordinates {r,9). This is unlike the standard bumpy 
formalism that introduces six non-vanishing metric com- 
ponents that depend on four arbitrary functions (cour- 
tesy of the Newman- Janis algorithm applied to a non- 
Kerr metric). Furthermore, this analysis shows that two 
of the six arbitrary functions introduced in the general- 
ized bumpy scheme of the previous section must not be 
independent, as argued earlier. 

Following the insight from the previous section, 
we will henceforth allow these 5 metric components 
{h°'',h^^,h°^,h^^,h^^) to be arbitrary functions of 
(r, 9) , although technically there are only four indepen- 
dent degrees of freedom. We do this because it eases 
the analytic calculations to come when one investigates 
which conditions must satisfy for there to exist an 
approximate second-order Killing tensor. Moreover, it 
allows us to relax the undesirable requirement, implicit 
in Eq. ((54|) . that when SV ^ 0, then both htt and h^^ 
must be non-zero in the a — > limit. 



B. Existence Conditions for the Carter constant 



Let us now follow the same methodology of Sec. lIII Bj to 
determine what conditions the arbitrary functions must 
satisfy in order for there to be a second-order Killing ten- 
sor. We begin by parameterizing the Killing tensor just 
as in Eq. ([M)) . with the expansion of the null vectors as 
in Eq. (|25p and the replacement of BK DK every- 
where. With this at hand, the Killing equation acquires 
the same structure as Eq. ([^5]). 

We have solved the Killing equations with Maple and 
the GrTcnsorll package [s^ to find that the null vectors 
must satisfy: 



Sh" 



51" = 



A 



((S^oK + 71) +74,'5Zdk + 71, 



LDK 



('^^DK+7l)+73 



-^Sll, + ^,,511,, -^SF^^ + 73 



(55) 



where ji = ^i{r) are arbitrary functions of radius, while 
i5ZnK and h^f are arbitrary functions of both (r, 0). As 
before with Sk^^, we find below that the function 51^^^ 
does not enter the metric perturbation, so we can set it 
to zero. 

The functions 7^ are completely determined by the 
metric perturbation: 
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2a^r'^l\sm^0 p^{r'^ + a^)/^ a A sin^ 6* Pg^'^ 2A P°*^ 



„2 PDK ®3 ^2 pDK + ^2 pDK ^4 

p2A2 sin^ 6* d7i a A2(S + 2a2Mr sin^ 61) sin^ 61 ^73 A^ (p2 _ 4ii/r) sin^ 6^ ^74 



2M Pf"^ dr 2M Pf"^ dr 2M Pf""^ dr 



- -—I - —63 . 



(r2+a2)2 , A^ 2(r2 + a2)^,, , A^ 2A2 sin^ , 2A 



2 



= — + — '-hf^ + -e3 73 + — 74, 

a ^ a 

dh°^ 2r,„^ 2a2sin6lcos6l,„^ 2r 2r^ 

d'^K^ SaM sine cose P^"^ AaMr{r^ + a^) sine cose dh^^ 2a^ sin^ e P^"^ 2r Pf"" 

-Q^ - J8 re ji -p5^ 14, ~^pu^^t4> 

AaMrsin^eP^^ AaM r sin^ e P^"^ AaMr P^"^ 2sin^ 9 P^^ 

„4 DDK ^4 pDK^l ^ „4 pDK 3 + 4 pDK '^^ 

F ^16 F ^1 F ^1 r ^1 

IGaMsin^ePfi'^ 2a P^^ d-fi 2 sin^ 61 P^f ^73 2a sin^ 6* Pf/ ^74 aAsin^ 61^2^1 

Asin^g ,^ , ^ 2,,r • 2./'73 aA(p2 - 4Mr) sin^ ^^74 
4 — (S + 2a A/rsm 6')— ^ j — ^ , (56) 



where recall that p^ = — a?' cos^ 9, Qr = 0r(^) is an 
arbitrary function of polar angle, while P^"^ are polyno- 
mials in r and cosfl, given explicitly in Appendix [Xj No- 
tice that many of these relations are integro-differential, 
as I is defined as 



1. Non-Rotating Limit 

Let us now take the non-rotating limit, i.e. a — > 0, 
of the DK Carter conditions. As before, we expand all 
arbitrary functions as in Eqs. ((3T|) and ((32|) . Imposing 
the constraints 



I 



dr 



2a? sin 9 cos e 
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O , 2 0^71 

2r 71 + p — — 
dr 



(57) 



Notice also that the component h^f is free and thus there 
are here truly only four independent metric components. 

We see then that any metric perturbation that 
satisfies Eq. ([56)) will possess a Carter constant. If given 
a specific non-Kerr metric, one can then use these equa- 
tions to reconstruct the 7^ and 62 functions to automati- 
cally obtain the perturbative second-order Killing tensor 
associated with this spacetime. 

With this at hand, one can now construct the constants 
of the motion. The metric in Eqs. (|53p and (j54p remains 
stationary and axisymmetric since all metric components 
depend only on (r, e). Therefore, the constants of the 
motion can be expanded in exactly the same way as in 
Eqs. ([55]) - ([55)1 . With this at hand, the geodesic equations 
can be rewritten in first-order form in exactly the same 
way as in the BK case, except that here different metric 
components are non- vanishing. Of course, all perturbed 
quantities must here have a DK superscript. 



63,0 

73,0 
74,0 



0, 
0, 



LDK 

'htfl 



71,0 



73,1 



LDK 

"■rr,0 



LDK 

"■rr,0 
~2^ 



of 



2r3/ 



4M^oK 
2 tt,o 



2/2 2 ' 
forces the metric components to take the form 



. DK 



. DK 
''tt,0 I 



. DK 



(58) 



0. (59) 



All other pieces of 0{a'^) in the remaining metric com- 
ponents can be set to zero by setting {h°g, hgg, hf^) to 
zero. With these choices, the differential conditions are 
automatically satisfied, where we have set all integration 
constants to zero. 

The Killing tensor is then given by Eq. ((24| . with the 
perturbed null vector components 



^^'dk 



51" = 

"'-DK 



"rr,0 



LDK 

"■rr,0 



+ 



"■tt.o 



2/2 

2 ,0,0,O(a) 



(60) 



LDK 
"tt,0 



2/ 
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where here we have set 51^^ = 0. As before, the Kilhng 
tensor reduces exactly to that of the Schwarzschild 
spacetime, with the only non-zero components being 
(60,^00) =r-4(l,sin2 0). 



RELATING PARAMETERIZATIONS 



To Each Other 



In the previous sections, we proposed two different pa- 
rameterizations of deformed spacetimes suitable for mod- 
eling alternative theory predictions. These paramcteri- 
zations are related via a gauge transformation. Under a 
general diffeomorphism, the metric transforms according 
to h'j^^ — = + V(^^y), where we parameterize 
the generating vectors as 



(61) 



The question is whether a generating vector exists that 
could take a generic metric perturbation to one of 
h^'^ form. The DK parameterization has an (r, 6) com- 
ponent that is absent in the BK one, while the BK one 
has {t, r) and (r, </>) components that are absent in the 
DK parameterization. 

Let us assume that we have a certain metric in DK 
form. The first task is to remove the (r, 9) component, 
i.e. to find a diffeomorphism such that = 0. This 
can be achieved by requiring that 



1^ 
A de 



96 
dr 



= 0. 



(62) 



whose solution is 



Ci=Fi(r)-A 



96 
dr 



(63) 



where Fi{r) is a free integration function and ^2(1", I 
free. 



The generating vector of Eq. (|6T|) with the condition 
in Eq. not only guarantees that h'^g = 0, but 

also ensures that the only non- vanishing components of 
the gauge-transformed metrics are the (i, t), {t, r), {t, (p), 
(r, r), {r,(p), {9,9) and (0,0), exactly the same non-zero 
components as in the BK metric. The new components 



are 



7 DK 

'Hr 



, DK 



"-rr 



,OK Mp^ 



2a? Mr sin 61 cos 6* 



-6, 



1/ 2Mr\dio aMr sin^ 9 d^3 
•^td, + ?i 



S2 

2aMr{r^ + a^) cos 9 sin 9 



S2 



6, 



a r a smt)cost) 
A 



A dr 



2S 



2aMr 



dr 



+ a^Asin^^l^l 
dr J 



= hgg+r^l- 0^81119 0089^2 + ^^, 



89 



J DK 



cos^ 9 [2r (r^ + M^) + AM] 
a'^ cos'^ 9{r - M)} £_i + sin 6* cos 



A- 



2A/r(r2 -|- , 



6, 



(64) 



where recall that = r^ — cos^ 9 and 6 is given by 
Eq. dMl)- 

The above result can be simplified somewhat by setting 
6 = and Fi{r) = in Eq. ([63]). Notice that the (i,r) 
and (r, 0) components are not modified by these vector 
components. The modified components then become 



"■tt 



LDK 
"t0 



ldk' 
ldk' 



LDK 
"tt 



LDK 
"t0 



LDK 



Mp" 



a, 



S2 

aMp^ s\-D^9 



a?r 



r6 



sm 



6 

A dr 



c? cos 



cos 



E2 
2 9 [2r (r^ 

4 



M)}6- 



AAf 



(65) 



We have thus found a generic diffeomorphism that 
maps a DK metric to one that has the same non-zero 
components as a BK metric. Notice that 6 S'Hd 6 o^ily 
enter to generate (t, r) and (r, 0) components. If we know 
the form of the components that we are trying to map to, 
then we could solve for these two vector components. For 
example, let us assume that h\l^ and are given and 

we wish to find 6,3 such that h^^^ = hf^ and h^^' = h^^. 
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This implies 



, 2aMr , 
2S dr 



dr 



(67) 



Let us, for one moment, allow ^0,3 to be arbitrary func- 
tions of (r, 6) . The solution to the above system is then 

[/if," cos2 9 {2aHlr - - a*) 
- hf^ {2a^Mr + a^r^ + a^) - 2Mrah^f\ , (68) 



2aAfrsin^ 6'/if^'' 



(69) 



But of course, to avoid introducing other spurious com- 
ponents to the DK metric, such as hf^, we must require 
that ^0.3 be functions of r only. This implies that the 
integrands must be themselves also functions of r only. 
Setting the integrands equal to F2 (r) and F3 (r) , we find 
the conditions 



'Hr ~ 



7 BK _ 



F2{r) {p^ - 2Mr) + 2Mras\n^ OFj,{r) 



2Mra sin^ 0F2{r) 



y + 2Mr)] sin2 6'. 



(70) 



Provided hf^ and h^^ components can be written in the 
above form, the generating function becomes 



F2ir)dr, 



F3{r)dr. (71) 



B. To Alternative Theories 

1. Dynamical Chern- Simons Modified Gravity 

The BK and DK parameterizations can also be mapped 
to known alternative theory BH metrics. In dynamical 
CS gravity [l3|, one can employ the slow- rotation ap- 
proximation, where one assumes the BH's spin angular 
momentum is small, l^l/M^ ^ 1, and the small-coupling 
approximation, where one assumes the theory's correc- 
tions are small deformations away from GR to find an 
analytical BH solution. Yunes and Prctorius [SQj found 
that this solution is simply g^^, = gfj^^ + K^^, where g^y 
is the Kerr metric (expanded in the slow- rotation limit), 
while the only non- vanishing component of the deforma- 
tion tensor in Boyer-Lindquist coordinates is 



LCS 
%0 



5, a . 2 

— Crq 1- Sm 



1 



12 A/ 27A/^ 



7r 



10?" 



(72) 



where a is the Kerr spin parameter, M is the BH mass 
and is a dimensionless constant that depends on the 
CS couplings and is assumed to be small. 

Such a non-GR BH solution can be mapped to the gen- 
eralized bumpy parameterization of Eqs. (|2ip and (j30p 
via §2 = Sg^^ = 65, and all other Si and Oi vanish. The 
quantity Sg'^^ is given by [4lj 

^ , .70r2 + 120Afr+189Af2 
= -CosaM^ n2.B(l„2A/A) ' ^''^ 

With these choices, the full null vectors become 

r = 



.-1,0, — + Sg^ 



(74) 



which agrees with Eq. (35) of [4l|. Moreover, with these 
choices all the metric components vanish except for the 
{t, 0) one, which must be equal to Eq. ([7^ for the differ- 
ential constraint in Eq. (jSOj to be satisfied. 

Similarly, we can also map it to the DK parameteriza- 
tion of Eqs. ([55)1 and ([5^ via 73 = Sg^^, and all other ar- 
bitrary functions vanish. With this, the perturbation to 
the principal null directions are equal to Eq. ((74)) . which 
agrees with [4l|. Moreover, we have checked that Eq. ([7^ 
satisfies the differential constraint found in Sec. IIVI 



2. Dynamical, Quadratic Gravity 

We can similarly map the BK and DK parameteriza- 
tions to the BH solution found in dynamical quadratic 
gravity [25| . Treating GR corrections as small deforma- 
tions (i.e. working in the small-coupling limit), Yunes and 
Stein [25! found that the unique non-spinning BH solu- 
tion is (7^1/ = ^^1/ -I- /iJJ^, where g^^^ is the Schwarzschild 
metric, while h';^^ in Schwarzschild coordinates is given 
by M 



'ht 



Cqg Af3 



SOAP 



1 



26M _^ 66 AP ^ 96 



_CqgM2 

j-2 y,2 

368 Af ^ 



M 

r 



52 M 



r2 



2AP 



16 Af4 



(75) 



where Ad is the mass of the Schwarzschild BH and Cqg is 
a dimensionless constant that depends on the couplings 
of the theory and is assumed to be small 

For a spherically symmetry background, one can easily 
show that the Schwarzschild Killing tensor (i.e. one whose 
only non-vanishing components are S,ee = and = 
r"* sin^ 9) satisfies the Killing equation, regardless of the 
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functional form of the (t, t) and (r, r) components of the 
metric. This means that the perturbed vectors Sk" and 
SI" must adjust so that the above is true. We find that 
the following vector components do just that: 



5¥ = 



51' 



2/r6 
128 M3 
5 

16 368 A/5 

r2^4 



8 A/ 14A/2 
6 r 



7\/ 52 A'/2 2Af3 
1 + — + ^ + — ^ 



2/2r6 
48 448 Af 5 



4 AT 26Af 2 32 
1 + -— + 



(76) 



and all other components vanish. 

That perturbed null vectors exist to reduce the Killing 
tensor to its Schwarzschild counterpart must surely be 
the case, as we have already shown in Sec. HUB II 
and IIVB II that in the non-rotating limit both the BK 
and DK parameterizations allow for generic (t, t) and 
(r, r) deformations. 



3. Arbitrary perturbations: perturbation to 73, 
perturbations to 71 and 74 

We can understand the functions that parametrize the 
metric perturbation by considering perturbations where 
we allow only a few of the parametrizing functions to 
be nonzero. Consider a perturbation in the deformed 
Kerr parametrization of Eqs. ((53|) and ((54|) described by 
7i = 74 = O3 = and /i^e = hgg = 0. In the small a 
limit, the only nonzero component of the metric is the 
(t,(f)) component: 

ht4, = -f{r^smH^3{r) + C), (77) 

where C is a constant. Now consider a perturbation in 
the DK parametrization whose only nonzero parameters 
are 71 and 74. In the small a limit, this produces pertur- 
bations to the metric given by 

/7i(r) + 2/274(7-), 
4Af)d74" 



htt 



f sin^ e 



2r274(r) 



2A/ 



dr 



(78) 



We then see that in the DK parameterization 73 con- 
trols modifications to frame-dragging (and thus the er- 
gospherc), while 71 and 74 control modifications to the 
location of the event-horizon and the innermost-stable 
circular orbit. 



VI. DISCUSSION 

The detection of GWs from EMRIs should allow for the 
detailed mapping of the spacetime metric on which the 
SCOs move. Such a mapping in turn allows for null tests 
of GR, as one would in principle be able to constrain 
deviations of the background spacetime from the Kerr 
metric. To carry out such tests, however, one requires 
a parameterization of the metric tensor that can handle 
model- independent, non-GR deviations. This was the 
purpose of this paper. 

We have constructed two such parameterizations: one 
as a generalization of the bumpy BH formalism (the BK 
scheme) and one as a generalized deformation of the Kerr 
metric (the DK scheme). The former promotes the non- 
vanishing components of the metric perturbation that 
contain bumps to arbitrary functions of radius and polar 
angle. The latter takes the most general stationary, ax- 
isymmetric line element without assuming Ricci-flatness 
and constructs a generic metric deformation from it. In 
both cases, the arbitrary functions introduced are con- 
strained by requiring that there exists an approximate 
second-order Killing tensor. 

These schemes differ from each other in the met- 
ric components that are assumed not to vanish. We 
have found a gauge transformation, however, that relates 
them. That is, we have shown that a generating vector 
exists such that the BK metric can be mapped to the 
DK one. We have also mapped both parameterizations 
to known analytical non-GR solutions, thus automati- 
cally finding their respective Killing tensors and Carter 
constants. 

The perturbative nature of our approach puts some 
limitation on the generality of the deformation. All 
throughout we restricted the metric deformation to be 
a small deviation away from the GR solution. But due 
to the structure of the perturbation, there can be re- 
gions in spacetime where this deformation dominates over 
the GR background. For example, the full metric tensor 
might lose Lorentzian signature if the metric perturba- 
tion dominates over the Kerr background and it possesses 
a certain definite sign, i.e., the metric will no longer have 
det(g) < 0. Clearly, such regions are unphysical, and 
the coupling constants that control the magnitude of the 
perturbation should be adjusted appropriately to ensure 
that if they exist, they are hidden inside the horizon. 

We should emphasize that we were guided by two prin- 
ciples in the construction of the parameterized schemes 
proposed here. First, we wanted to ensure that the pa- 
rameterizations would easily map to known analytic so- 
lutions. Second, we wished to retain a smooth Kerr limit, 
such that when the deformation parameters are taken to 
zero, the deformed metric goes smoothly to Kerr. This in 
turn guaranteed that certain properties of the Kerr back- 
ground were retained, such as the existence of a horizon. 
We additionally required that the metric tensors allowed 
a second-order perturbative Killing tensor, so that there 
exists a perturbative Carter constant. This latter re- 
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quircmcnt is not strictly necessary to perform null tests, 
as one could build GWs from the evolution of the full 
second-order equations of motion without rewriting them 
in first-order form. 

An interesting avenue of future research would be to in- 
vestigate the Petrov type [42| of the BK and DK metrics. 
Brink has shown that GR spacetimes that admit a Carter 
constant must also be of Type D . In alternative the- 
ories of gravity, however, a formal mathematical proof of 
the previous statement is lacking. We have here found 
Ricci-curved metrics that admit a Carter constant. One 
could now attempt to prove that the perturbative Carter 
existence conditions found in this paper also automati- 
cally imply that the spacetime is approximately Petrov 
Type D. This is a purely mathematical result that is be- 
yond the scope of this paper. 

An interesting physical question that could now be an- 
swered with either the DK or BK frameworks is whether 
LISA could place interesting constraints on non-GR Kerr 
deformations. To answer this question, one would have 
to evolve geodesies in either framework and construct 



the associated GWs. One could then carry out a Fisher 
analysis to determine the magnitude of the non-GR de- 
viation from Kerr that could be measured or constrained 
by LISA. We will investigate this problem in a forthcom- 
ing publication. 
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Appendix A: Polynomial Functions 



We here provide explicit expressions for the polynomi- 
als appearing in the conditions that guarantee the exis- 
tence of a Carter constant. In Eq. ((30|). the polynomials 
Pf"^ are given by 
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Pf = + ba^r^ + 2aV^ - 2a^ cos^ e{2r^ + a^r'^ + a^) - cos^ 9{r^ - a^) , 

PJ"' = r^M + c?r^{r + 2M) - cos^ B{r'^{r + 2Af) - a2(r - 2A/)) - a" co%^ e{r - M) , 

P^"^ = r{r^ +a^{r + 4:M)) + a^cos^9{r^ -Ulr + a^), 

P^"^ = r^{r^{r - 2M) + a^{3r - 4M)) - 2a^r^ cos^ 9{r^ - 2Mr - a^) - cos'' Qi^^r^ - 2Mr + c?) , 
Pg^K = (r - 2M){r^{r + 2M) +a'^{r + AM)) - cos^ ^(r^ - AM'^r'^ - a^{r'^ - 6Mr + 16M^) - 2a'') 

-a^ cos'' e{2r^{r^ - 2A/r + AM^) + a^r{r^ + AMr - AM^) - a'^{r - 2M)) - cos^ e{r - 2M)A , 

3^ lA{p^A-2Mrp^+AMra^+AMr^)sm^e I a A {sin (9))^ {-p^ + AMr) 

P, ="2 7 ' ^7 =2 ^3 ' (Al, 

Pg^^ - (8 r" - 4 r^p" + 16 r^A/^ + 2 p^ - 44 a^Mp^A - 8 a'°M - 68 a^Mr^p^ ^^28 a^M'^r^p^ - 56 a^rM^p^ 

- 12 a^r^p^ - 16 r^a^p^ - 6 a^rp^ _ 36 a^Mr^p^ + 40 a^Mr^p^ + 32 r^a^ - 3 r^p^ - 8 ri°Af - 20 o^rM^p^ 

- 16 r^Ma^ - 16 r^'a^M^ - 48 r^a^M^ + 16 r^Ma^ - 56 a^M'^r^p'' + 64 a^M^r^p^ + 72 a^rM^p^ + 43 r^a^ 

+ 32 r^a^ + 8 r^a^ + 60 a^Mp^A - 18 a^Mp^A - 80 a^M^r^ + 8 a® A/A + 16 a^M'^r + 36 a^Mp^ - 32 r^A'/2p2 

- 4 r®Afp2 - 12 r^a^p^ - 12 - 4 Mr^ p^ - A a^rp* - AA a^Mp^ + 12 Mr*p^ - 8 M^r^p^ + 24 M^r^p^ 
+ 2r^a'^p^ + 5 a^rp^ + 12 a^Mp^)sm'^ 9, 

Pg^« = (4 r*^ + 2 p^a^ - 4 a^p^ A - 8 a^p^Mr + 2 a^p^ + 8 r^c? - p" A - 2 p^Af r - 2 r^p^ + 4 a^A + 8 Mra" - A a^) p" , 

P^^, = a sin^ 9 (-16 A/^V^ + 8 A/V^p^ - 4 A/r^p" - 20 Mr^ - 32 a^A/V^ + 24 r^M - 24 a*r^ - 8 a® A// - 8 

+ 56 a^Mr^ + 40 a^Af r" - 24 r^a^ - 8 r^a^ + 8 a^A/A + 16 a^M'^r - 12 a^p^M + 28 a^Afp^ - 36 a^Afp^ A 

+ 12 aV^Af A - 48 a^M^r^ + 8 ra'^M^p'^ - 3 a\p^ + 16 r^A/^ + 2 a^A/p*^ + 18 a^r^p^ + 6 r^p^ + 18 r^a^p^ 

+ 6 a^rp^ + + 48 a^M^r^p^ ~ 40 a^M^rp^ - 48 a^Afr^p^) , 

Pg^^K ^ p^[Ar^ + 2 p'^c? - A a^p^A - 8 a^p^Mr + 2 a^p^ + 8 r^a^ - p-* A - 2 p^A/r - 2 r^p^ + 4 a-'A + 8 A/ra"* - 4 a^) , 

PfJ^„ = sin^ (-12 rA/a^p'' - 48 a^r^ + 6 r^a^p^ - 48 a^r^ + 16 a'^Af ^A + 32 a^M^r + 8 p^ A - 8 a^Af ^p^ - 16 r^^ 

- 16 A/V^ + 6 a^p^Mr + 24 r^a^p^ + 32 A/ - 32 a^p'^Mr^ - 8 + 8 a^M'^p^ + 40 M^r^ p^ + 16 rMa^p^ 

- 32 r^Ma^p^ + 8 r^p^ - 16 a^Af^ + 96 a^AfVp^ - 16 r'^a^ + 32 r'^Afa" - 80 r^a^Af^ - 48 r^a^Af ^ + 16 r^Ma^p^ 
+ 24 aV'^p^ - 6 Mr^p^ + 3 r^p*^ + 2 r^p^ + 5 a^p^ + 64 r'^A/a^ + 24 Af r^" - 24 AfV^p^ - 48 Mr'^ p^ 

- 24a^Af2p2sm2 6i-6aV*^A + 6aV^A-4aV) , 

Pio,d = P** (4 + 2 p^a^ - 4 a^p^ A - 8 a^p^Mr + 2 a^p^ + 8 r^a^ - p** A - 2 p^Afr - 2 r^p^ + 4 a^A + 8 Mra^ - Aa^) , 

Pfj^^ -8 aM sin^ 6* (r^p" - 5 pV^ + 4 r'^ - 3 ra^p^ + 6 aV'V^ + aV''A/ + a^p^MA + 2 a^p^M^r - 3 a^Mp^ - 8 a^Mr* 

+ 3 a*?y - 2 A/r*^ - 4 a^r^ + 3 A/r^p^ + 2 a''A/A + 4 a^AfV - 2 a^Af) , 

Pfj^^^ = p4 (4 + 2 p^a^ - 4 A - 8 a'^p^Mr + 2 a^p^ + 8 r^a^ -p^A-2 p^Mr - 2 r^p^ + 4 a^A + 8 A/ra^ - 4 a^) , 

ATri = -2 r (-28 a^p^ A - 56 a^p^Mr + 18 a*p^ + p^ + 2 p** A + 4 p''A/r + 4 p^a^ - 18 r^p^ + 16 

+ 32 r^a^ + 16 a* A + 32 Af ra^ - 16 a^) , 

Pi2,d = p^ (4 + 2 p^a^ - A a^p^A - 8 a^p^Mr + 2 a^p^ + 8 r^a^ - p" A - 2 p^Afr - 2 r^p^ + 4 a^A + 8 Af ra^ - 4 a^) , 

Pfg^^^ = 2 (sin (61))^ (8 a^p^A + 16 aV^'^^ - 8 aV -P^- 6 p^A- 12 p''A4"r + 8 p''a2 + 24 rV- 16 16 r''a2) , 

Piz.d = p^ (4 + 2 p^a^ - 4 a^p^ A - 8 a^p^Mr + 2 a^p^ + 8 r^a^ - p" A - 2 p'^Mr - 2 r'^p'^ + 4 a^A + 8 A/ra^ - Aa^) , 

_ -4 aAfr (sin (61))^ (-4 r'' + p^A + 2 p^Afr- 4 a2p2 + 4 r^a^) 

" ^ p4 (4 r6 + 2 p-^a^ - A o?p^A - 8 a^p'^Mr + 2 a4p2 + 8 r'^a? - p^A - 2 p^Mr - 2 r*p'^ + A a^A + 8 Af ra^ -Aa^)' 

I 



Similarly, the polynomial functions that appear in the 
Carter conditions for the DK metric of Eq. ((55)) are 
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P^^ = r^{r'^ + + + 2a'^) - 2a^ cos^ e{2r^ + + a^) - cos^ e{r^ - a^) , 
^2°" = r^ir^AI + a^r + 2a^M) - cos^ e{r'^{r + 2M) -a^{r- 2M)) - a" cos" d{r - M) , 
= r^{r - 2M){r'^{r + 2M) + a^{r + 4M)) + 2r^a^ cos^ Oir"^ - 2Mr + AM"^ + a^) 
W COS-* 61 (r - 2Af )A , 

pDK ^ j,3(^2(^^ _ 2A/) + a2(3r - 4A/)) - 2r'^a^ cos^ 6l(r2 - 2Mr - a^) - a"^ cos'' e{Zr^ - 2Mr + a^) , 
^5°'" = r'^(r^ + ia^r^ + SaV^ + 2a^) - Sa^r^ cos^ 6l(3r2 - a^) + a"* cos'' 6l(5r^ - SaV"* + ea'^r^ + 2a^) 

cos^ 6l(2r4 - Sa^r^ - a'') , 
Pg^-^ = r4(r2 - 6a2) + Sa^r^ cos^ 6'(3r2 + Aa^) - a" cos^ 6'(9r2 - 2a^) - a*^ cos^ 9 , 
P^"^ = + lOaV'' + Qa^r"^ - cos^ 6'(llr'' + IBaV^ + 10a'') + a" cos^^ 6'(5r2 + 60^) + cos^ 6* , 

Pg^-^ = r2(3r4 + 5aV2 - 2a^) ~ 2c? cos^ ^(r^ + c?)[2r^ - Sa^) + cos^ - Sa^) , 
P^^ = r2(3r2 - a^) _ 3 cos^ 6l(r'' - a'') + cos" 61 (r^ - Sa^) , 

Pfi"" = -r''(3r6 - 2Mr5 + 240^"* - ISa^Mr^ - Sa^M^^ + lOa^r^ - 24a''A/r + A&a^M'^ + 6a^) 

W cos^ 6'(21r^ - 18A//r^ - 8MV^ + 18aV^ - \2c?Mr' + 72a2A/V + SSfl-^r" - 'i2a^Mr' + 24a''A//V2 
+12a^r2 - lOa^Af 2) + a" cos** ^(llr^ + %Mr^ - 24:A4\^ - Vla^r^ + 22c?Mr^ - 2Ac?M'^r^ + Za^"^ 
-2Aa^Mr + 240* A/^ + 2a^) + cos^ 6l(3r'' - QMr^ - OaV^ + ISa^A/r - Sa^Af ^ _ a'^) , 

AT = r^{r^M + SaV^ - 6a2A//r - a") - cos^ 6l(3r5 - 3A/r'' - Sa^Afr^ - 'ia^r + 2a^M) 
cos" e{r^ - 3aV + a^M) , 

P""^ r3(r'* + 9aV"-10a2Afr3 + 6a"r2-8a''Afr + 2a^) 

- cos^ er{r^ + ISaV^ - 22a^A.fr^ + 15a"r" - 26a" A/r^ + 20aS'2 - 120*^ Afr + 6a*^) 
W cos" 6l(9r^ - 12Afr^ + 6a2r-^ - 22a2Af r" + 27a"r3 - 18a"A/r2 + Qa^r + 40*^ Af) 
+a" cos^ 6l(3r''^ + 4A-f r" - lOaV^ + Qa^Mr"^ + 3a"r - 6a" Af) + aP cos* 0{r^ - Sa^r + 20^ M) , 

Pit = ^^(37-*^ - 4Afr^ + 22aV^ - 30a2Afr^ + 8a^AfV" + 29a"r" - 14a" Afr^ + 12a'^r^ - 120*^ Afr - 16a^A/2 + 2a*) 
-a^ cos^ 6l(15r'' - 24A//r* + 8A/V^ + lOa^Afr^ + 21a"r5 - 34a" Afr" - 32a" AfV^ + 24a^r3 + 12a^Afr2 
-32a^A//V + 6a*r + 8a*Af) - a" cos" 6l(13r^ - 32Afr^ + 6a^r^ + Ua^A/r" + ?>a'^M'^r^ - ia^r^ - 'i2a*Mr'^ 
+48a"AfV + 4a^r - lOa^A/) - cos^ 6i(3r^ - 12Afr" + 8A/2r3 - 2a^r^ + ISa^Mr^ - IQa^M'^r - 5a"r + 6a" Af) , 

P°^ = r^{r^ - M'l'^r^ + 9a2r" - 22a^Mr^ + 6a"r2 - 4a"A/r + 16a" A/^ + 2a^) 

-a^ cos^ 6l(9r^ - 28A'/r^ + 9aV5 - 18a2Af r" + 16a2Af + lea^A//^^ + 18a"r3 - ma^Mr^ + 6a^r - 8a^Af ) 
-a" cos" 6lr(3r" - 4A//r3 + (8Af ^ - 9a2) + \%a^Mr - \%c?M'^) - cos^ e{r^ - 3a^r + 2a^M) . (A2) 



Appendix B: Equations of Motion in Alternative 
Theories 

The derivation of the equations of motion in Sec. Ill Al 
is fairly general. In particular, this derivation is indepen- 
dent of the metric used; in no place did we use that the 
spacetime was Kerr. We required the divergence of the 
test particle stress-energy tensor vanishes; this condition 
is always true in GR, because of a combination of local 
stress-energy conservation and the equivalence principle. 

One might wonder whether the divergence-free condi- 
tion of the stress-energy tensor holds in more general the- 
ories. In any metric GR deformation, the field equations 
will take the form 

Gap + 'HaP = T^p + T^^p , (Bl) 
where Gap is the Einstein tensor, "Hap is a tcnsorial de- 



I 

formation of the Einstein equations, and is a possi- 
ble stress-energy modification, associated with additional 
fields. The divergence of this equation then leads to 

V'^Uap - V"C^-' + V"T^^ , (B2) 

since the Bianchi identities force the divergence of the 
Einstein tensor to vanish. The Bianchi identities hold in 
alternative theories, as this is a geometric constraint and 
not one that derives from the action. We see then that 
the divergence of the matter stress-energy tensor vanishes 
independently provided 

V'^Hap = V"r«! • (B3) 

Whether this condition [Eq. (jB3[) ] is satisfied depends 
somewhat on the theory of interest. Theories that in- 
clude additional degrees of freedom that couple both to 
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the geometry and have their own dynamics usually sat- 
isfy Eq. (|B3p . This is because additional equations of 
motion arise upon variation of the action with respect 
to these additional degrees of freedom, and these addi- 
tional equations reduce to Eq. (|B3p . Such is the case, for 
example, in dynamical CS modified gravity [4l|- If no 
additional degrees of freedom arc introduced, the satis- 
faction of Eq. (|B3|) depends on whether the divergence of 
the new tensor T-Lap vanishes, which need not in general 
be the case. 

Recently, it was shown that the equations of motion 
are geodesic to leading order in the mass-ratio for any 
classical field theory that satisfies the following con- 
straints [l6j : 



• It derives from a diffcomorphism-covariant La- 
grangian, ensuring a Bianchi identity; 

• It leads to second-order field equations. 

The second condition seems somewhat too stringent, as 
we know of examples where third-order field equations 
still lead to geodesic motion, i.e. dynamical CS grav- 
ity [4l|. Therefore, it seems reasonable to assume that 
this condition could be relaxed in the future. Based on 
this, we take the viewpoint that the equations of motion 
are geodesic even in the class of alternative theories we 
consider here. 



